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Abstract 



The production fractions of charged and neutral 6-hadrons in 6-quark events 
from decays have been measured with the DELPHI detector at LEP. An 
algorithm has been developed, based on a neural network, to estimate the 
charge of the weakly-decaying 6-hadron by distinguishing its decay products 
from particles produced at the primary vertex. From the data taken in the 
years 1994 and 1995, the fraction of 6-quarks fragmenting into positively charged 
weakly-decaying 6-hadrons has been measured to be: 

/+ = (42.09 ± 0.82(stat.) ± 0.89(syst.))%. 

Subtracting the rates for charged and baryons gives the production 
fraction of mesons: 

/^^ = (40.99 ± 0.82(stat.) ± l.ll(syst.))%. 



(Accepted by Phys. Lett. B) 
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1 Introduction 

The branching fractions of the 6-quark into the different species of &-hadrons are an 
important input and source of systematic uncertainty for many measurements in the 
heavy flavour sector where 6-hadrons are produced in jets, e.g. analyses on S-meson 
oscillations or CKM elements at LEP. Furthermore, these production fractions give insight 
into the fragmentation process. Since 6-quarks at LEP are mainly produced directly in 
the decay of the boson, with negligible contributions from later processes like gluon 
splitting g bb, 6-hadron production fractions are sensitive to a certain step in the 
fragmentation process, namely the beginning of the fragmentation chain. This is not the 
case for analyses investigating inclusive particle production rates of hadrons which do not 
contain a primary heavy quark. 

The 6-hadron production fractions are defined as the probability of a b- or 6-quark 
to fragment into the corresponding 6-hadron: /^^ = BR{b B^) = BR{b — > 

fs, = BR{b ^ fiO) = BR{b ^ fB, = BR{b ^ B^) = BR{b ^ fi^), h-Mryon = 
BR{b — > anti-6-baryon) = BR{b — > 6-baryon). Furthermore, the production fractions for 
charged and neutral 6-hadrons are defined as = BR(^ X^) — BR(b — > X^) and 
f = BR{b X^) = BR{b X%), where X^, X^ and stand for any positively 
charged, negatively charged or neutral 6-hadron, respectively. With these definitions, 
fsi is also the production fraction of the 6-hadron type Bi, particle or antiparticle, in 
66-events. 

A direct measurement of these production fractions using exclusive decays is difficult, 
since there are many decay channels with small branching fractions having large relative 
uncertainties [1]. For the determination of /^j, , the inputs used are measurements of the 
product branching ratio BR{b — > 5^) ■ BR{B'^ D^l^uX) at LEP [2], measurements 
of the ratio fsJifsu + fsd) using events with exclusively reconstructed charm particles 
in semileptonic ^-decays or double semileptonic decays from CDF [3] and the mixing 
parameters x and Xd- The integrated mixing probability x, in an unbiased sample of 
neutral i?-mesons, has contributions from B^- and i?°-mesons^: % = fs^Xd + fs^Xs, 
where Xd and Xs are the integrated mixing probabilities for and S°-mesons. This 
allows the extraction of fs^ quite precisely. The baryon rate is estimated from similar 
product branching ratios, using A^/~ and correlations [3,4], and a measurement of 

proton production in 6-hadron decays [5]. No direct measurements of Jb^ or fs^ have 
been published so far. The averages for weakly-decaying 6-hadrons are hsted in [1] . The 
following assumptions are made: the 6-hadron production fractions are the same in Z — > 
bb decays at LEP and in high-p^ jets at the TEVATRON, /b, + fb-baryon + fsu + fsa ^ i 
and fs^ = Jb^- The latter two are applied as constraints in the averaging procedure. 
The combined result is /b, = (10.6 ± 1.3)%, h-baryon = (11-8 ± 2.0)% and Jb, = Jb^ = 
(38.8 ± 1.3)%. 

In Figure 1 a schematic picture of the 6-hadron production process is shown. One 

has to distinguish between the fractions of &-hadrons which are directly produced in the 
non-perturbative fragmentation process, denoted // in the following, and the fractions of 
weakly-decaying 6-hadrons, /j. Strong decays of primarily produced 6-hadrons can lead 
to // ^ /i, e.g. the presence of orbitally excited S**-mesons with their expected decays 
BT B^^U\ liavo the consecjuonce /^^ > /b^, /^^ < /b„ and /^^ < /b^. However, the 

^ Since x is mainly measured using leptons, the rates fs^ and /s^ have to be weighted by ratios of lifetimes, Tb^I^b 
and tb^ /tb , respectively. This has been omitted in the formulae for simplicity. 
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Figure 1: Schematic picture of the production mechanism of 6-hadrons. The rates, given 
in percent, of hadrons primarily produced in the fragmentation are denoted by /', the 
ones which decay through weak interaction (indicated by dashed hues) by /. Strong and 
electromagnetic decays are indicated by solid and dotted arrows, respectively, together 
with their (expected) branching ratios (for strong decays only single pion and kaon tran- 
sitions have been considered). The given rates are taken from simulation (Jetset 7.3 
Monte-Carlo model with parton shower option and parameter settings according to the 
DELPHI- Tuning [7]). The parameters giving the suppression of ss pairs and diquarks 
are, respectively, 7^ = P{ss)/P{uu) = P{ss)/ P{dd) =0.28 and P{qq)/ P{q) =0.1. 



equality /^^ = /^^, which originates from isospin symmetry, remains valid for weakly- 
decaying 6-hadrons (/b„ = /bJ.^ 

For the rates of charged and neutral 6-hadrons, an efficient algorithm has been de- 
veloped to distinguish charged particles from weak 5-decays, from their fragmentation 
counterparts produced at the primary event vertex. This allows an estimate of the charge 
of the weakly- decaying hadron to be made and thus a measurement of /"*" and fs^ 
can then be extracted from with small additional uncertainties. The data taken in 
the years 1994 and 1995, when the DELPHI detector was equipped with a double sided 
silicon vertex detector, have been used for the analysis. 

The simulation used the Jetset 7.3 model [6] with parton shower option and param- 
eters determined from earlier QCD studies [7], followed by a detailed detector simula- 
tion [8]. 

^In contrast to the D-system, the presence of S*-mesons does not change the rates of charged and neutral B-mesons, 

because their dominant decay mode is B* Bj. This is also the case for orbitally excited B**-mesons, if /g** = fB"", 

d ^ 

and isospin rules are used to calculate the dominant single pion transitions. 
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2 The DELPHI detector and event selection 

The DELPHI detector is described in detail in references [9,10]. The present analysis 
relies mainly on charged particles, measured using information provided by the central 
tracking detectors. 

• The micro Vertex Detector (YD) consists of three layers of silicon strip detectors 
at radii of 6.3, 9.0 and 10.9 cm. R(p coordinates^ in the plane perpendicular to the 
beam are measured in all three layers. The first and third layers also provide z 
information. The polar angle (9) coverage for a particle passing all three layers is 
from 44° to 136°. The single point precision has been estimated from real data to be 
about 8 fim in R(f) and (for charged particles crossing perpendicular to the module) 
about 9 fim in z. 

• The Inner Detector (ID) consists of an inner drift chamber with jet chamber 

geometry and 5 cylindrical layers of multiwire proportional chambers (MWPC). 
The jet chamber, between 12 and 23 cm in R and 23° and 157° in 9, consists of 24 
azimuthal sectors, each providing up to 24 R(j) points. From 1995 on, a longer ID 
has been operational, with polar angle coverage from 15° to 165° and replacing the 
MWPC by 5 layers of straw tube detectors. The precision on local track elements 
has been measured in muon pair events to be about 45 fim. in R(f). 

• The Time Projection Chamber (TPC) is the main tracking device of DELPHI. It 
provides up to 16 space points per particle trajectory for radii between 40 and 110 cm 
and polar angles between 39° and 141°. The precision on the track elements is about 
150 //m in Rcf) and about 600 fim in z. For particle identification a measurement of 
the specific energy loss (dE/dx) is provided by 192 sense wires located at the end 
caps of the drift volume. 

• The Outer Detector (OD) consists of 5 layers of drift tubes between radii of 197 
and 206 cm. Its polar angle coverage is from 42° to 138°. The OD provides 3 space 
points and 2 R(j) points per track. The single point precision is about 110 /xm in the 
R(f) plane and about 3.5 cm in the z direction. 

An event has been selected as multihadronic if the following requirements are satisfied: 

• There must be at least 5 charged particles in the event, each with momentum larger 
than 400 MeV/c and polar angle between 20° and 160°. 

• The total reconstructed energy of these charged particles has to exceed 12% of the 
centre-of-mass energy (assuming all particles to have the pion mass). 

• The total energy of the charged particles in each hemisphere (defined by the plane 
perpendicular to the beam axis) has to exceed 3% of the centre-of-mass energy. 

After these cuts, about 2 million events from the 1994 and 1995 runs have been re- 
tained. About 7 million simulated — > qq and 3.1 million — > bb events have been 
selected with the same cuts. 

Jets have been reconstructed with the LUCLUS algorithm [6] {djoin — 5 GeV/c) using 
charged and neutral particles. Two-jet events well within the acceptance of the vertex 
detector (| cos9thrust\ < 0.65) were selected. The event was divided in two hemispheres 
by the plane perpendicular to the thrust axis. 

The most important variables to tag or antitag bb events are the track impact param- 
eters of charged particles with respect to the primary vertex which is fit on an event-by- 
event basis using the position and size of the beamspot as constraints. From the track 

^In the standard DELPHI coordinate system, the z axis is along the electron direction, the x axis points towards the 
centre of LEP, and the y axis points upwards. The polar angle to the z axis is denoted by 0, and the aaimuthal angle 
around the z axis by 4>\ the radial coordinate is R = ^x'^ + y'^ . 
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impact parameters and their errors a probability is computed that a selected sample of 
charged particles originates from the primary vertex. To increase efficiency and purity, 
additional information e.g. from reconstructed secondary vertices and identified leptons 
are used. A combined discriminating variable is then used to select bb events. These 
methods are described in detail in [11]. The tagging of bb events was performed in the 
hemisphere opposite to the one which was used for the measurement. The cut on the dis- 
criminating variable of the combined 6-tagging has been chosen to give a 6- purity of about 
97.5%. A secondary vertex was fit in the hemisphere considered for the measurement. 
Hadronic interactions in the detector material were reconstructed using the algorithm 
described in [12]. Since the particle causing the interaction is lost in most of the cases, 
hemispheres with such interactions were rejected. 



3 Measurement of the rates of charged and neutral 
6-hadrons 
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Figure 2: The output of the neural network used to separate S-decay particles from frag- 
mentation particles. Shown are the data (closed circles), the simulation (solid histogram) 
and the contributions from weak 5-decay particles (dotted) and their fragmentation 
counterparts produced at the primary vertex (dashed). The latter distribution includes 
particles originating from strong decays of excited 6-hadrons. 



The basic idea for measuring the rates of charged and neutral 6-hadrons is to recon- 
struct the charge of the weakly-decaying hadron. Based on a neural network, for each 
charged particle in a hemisphere, a probability Pb that the particle originates from a 
6-hadron decay rather than from fragmentation is calculated. Charged particles are ac- 
cepted if their momentum exceeds 500 MeV/c and if at least one vertex detector hit has 
been associated. At least four charged particles had to be accepted in the hemisphere. 
The maximum number of charged particles in the hemisphere failing these acceptance 
cuts was limited to four. The input variables to the neural network are the probability 
that the charged particle track fits to the primary vertex, the momentum, the rapidity 
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Figure 3: (ctq^) (left) and cr{aQg) (right) versus the number of charged particles which 

have been used for the estimation of the vertex charge for data (closed circles) and 
simulation (histogram) before applying the correction as explained in the text. The 
ratios {aQg)data/ {<^QB)sim. and 0.1 ■ a {ag J data/ <^{<^QB)sim. are shown as dashed lines. The 
deviations are within ±3% over the whole range. 

of the particle with respect to the thrust axis, the reconstructed flight distance from the 
primary to the secondary vertex in the plane and its error. The last two are not spe- 
cific for the particles but for the hemisphere considered. They give additional information 

about the separation power of the other variables, especially the vertex probability. The 
distributions of the net output variable"^ are shown in Figure 2.^ 
A secondary- vertex charge is then constructed through 



where Nhem is the number of accepted particles in the hemisphere, Qi the charge of particle 
i and Pg,* its probability to stem from a 6-hadron decay as defined above. Assuming 
binomial statistics, an error on can be defined as 



This quantity does not account for particle losses due to inefficiencies in the track re- 
construction. (jQg is small if all charged particles are well classified, having values of Pb 
close to or 1, and gets larger the more particles have probabilities around 0.5. 

Parameters in the simulation possibly having an effect on the measurement have been 
adjusted to their measured values. They are discussed in detail in section 4 and fisted 
also in Table 1. 

*This variable can be interpreted as Bayesian a-posteriori probability. Pg is computed from this variable taking into 
account the ratio of b-hadron decay particles and fragmentation particles as taken from simulation. This is necessary 
because the neural network has been trained vi^ith equal numbers of charged particles from these two classes. 

^If not explicitely stated otherwise, all figures show the data from 1994 and 1995. 
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The shapes of the Qb distributions are directly affected by the number of charged 
particles which have been used in the reconstruction of the vertex charge. This number is 
larger in the data than in the simulation by 0.12 at a mean value of about 8.0. The shape 
of the particle multiplicity distribution is very similar in the data and the simulation. 
The simulated events are reweighted to get agreement in this distribution. The error of 
the vertex charge, defined in equation 2, also influences the Qb distributions because it 
is directly related to its width. Incorrect modelling of the shape of the Qb distribution in 
the simulation potentially biases the result. The dependence of the mean value and spread 
of o"Qg on the number of charged particles is shown in Figure 3. Data and simulation 
agree well, the deviations being within ±3%. The simulated events are reweighted to 
get these distributions in agreement, individually for any number of charged particles, to 
avoid any possible bias. To correct for the loss of charged particle tracks, the number 
of charged particles in the hemisphere not passing the track cuts is also brought into 
agreement. All these corrections are small. Comparing the Qb distributions of data 
and simulation shows slight shifts of the data distribution with respect to the simulation 
which is corrected depending on the absolute value of the polar angle 9 (these shifts are 
typically around 0.01 with the deviations from zero not being very significant). Such 
shifts can be caused by small differences in the material distribution of the detector in 
data and simulation causing a charge asymmetry. 




Figure 4: The fraction of 'opposite-sign' events versus Q^* fo^" 1994 (left) and 1995 
(right) comparing data (circles with error bars) and simulation (histogram), used for the 
calibration of Qb- For 1994, the solid (dashed) histogram shows the simulation before 
(after) applying the correction explained in the text. 



The vertex charge is sensitive to the charge of the 6-quark in the hemisphere, mainly 
in cases where charged 6-hadrons are produced. Since h- and 6-quarks are produced 
in pairs in Z decays, events where the vertex charge Qb can be determined in both 
hemispheres ('double tagged' events in the following) can be used to calibrate Qb on 
the data themselves. This is done in the following way. For 20116 events, where all 
cuts are passed in both hemispheres, 'opposite-sign' (called '05" in the following) and 
'same-sign' ('5*5") events arc defined through the vertex charges in both hemispheres: 
OS: Q^Q ■ Q\ < 0, SS: ■ Q% > 0- To be sensitive to the shape of the distributions, 
\Qb^\ > Qb* h^-s been required. The fractions of OS and SS events are directly related to 
the probability to tag the charge of the 6-quark in the hemisphere correctly. In Figure 4, 
the fraction of OS events, f^^, versus is shown for the 1994 and 1995 data sets. It can 
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be seen that the probabihty to tag the 6-quark charge correctly is shghtly overestimated 
in the simulation for 1994 (with a tendency to become worse when increasing Qs*, thus 
probing the wings of the distribution), whereas good agreement is found for 1995. The 
fraction of OS events is mainly determined by the Qb distributions of charged &-hadrons. 
To correct for the disagreement between data and simulation in 1994, the Qb distributions 
for positively and negatively charged 6-hadrons are reweighted in the simulation according 
to — 1 — (±ai • Qb — 0,2)^- The following parameters were found to make the fraction 
of OS events in the simulation consistent with the data over almost the full range of Q^* 
(see Figure 4): Oi = 0.31 and 0,2 = 0.34. Another, related, distribution is the fraction of 
double tagged events versus * with respect to all double tagged events. It has been 
verified that this distribution is in good agreement for both years of data taking (after 
applying the correction to the 1994 simulation). 




Figure 5: Distribution of the vertex charge, Qb, of the weakly-decaying b-hadron for the 
data (points with error bars) with the result of the fit superimposed (solid histogram) on 
a linear scale and on a logarithmic scale (as inlay). The shapes for neutral (dashed his- 
togram), negatively (dashed-dotted) and positively (dotted) charged 6-hadrons obtained 
from the simulation (in the fit, one single component was used for positively and nega- 
tively charged 6-hadrons) are also shown. The hatched histogram shows the contribution 
of non-66 events. 
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The measured Qb distribution has been fit by the corresponding shapes expected for 
charged and neutral 6-hadrons obtained from the simulation, while not separating the 
shapes for the positive and negative charges. 

A technique based on a binned log-likelihood method taking into account the limited 
statistics of the simulation has been used [13]. The non-66 background has been fixed 
to the value obtained from simulation. The real data distribution, corresponding to 
103.285 selected hemispheres, together with the fit result and the simulation prediction for 
neutral, positively and negatively charged 6-hadrons is shown in Figure 5. The result for 
/"'", the fraction of charged 6-hadrons in a sample of weakly-decaying 6-hadrons produced 
in Z° bb decays, is /+ = (41.84 ± 0.99(stat.))% for the 1994 data set and /+ = 
(42.65 ± 1.48(stat.))% for 1995 giving a combined value of 

/+ = (42.09 ± 0.82(stat.))%. (3) 

The result for f is given through /o = 1 - /+ = (57.91 ± 0.82(stat.))%, with /+ and f 
being fully anticorrelated. The P^r degree of freedom of the fits are 0.96 and 1.06 for 
the two years. 



4 Systematic checks and uncertainties 

Several cross-checks have been performed. The fit range and number of bins of the 
histograms used in the fit have been varied. The momentum cut has been varied in 
the range from 300 to 800 MeV/c and the maximum number of rejected tracks in the 
hemisphere between two and five. No significant change of the result has been found. 
The distributions of negatively and positively charged 6-hadrons have been fit separately. 
This gives 



5i?(6,6^X+) = (20.66±0.60(stat.jj7o 
BR{b,b^X^) = (21.16 ±0.59(stat.))%. 



for 1994 and 

BR(b,b^X+) = (21.37 ±0.87(stat.))% 
BR(b,b^Xs) = (21.29 ±0.86(stat.))%. 

for 1995. The two numbers are correlated (p^ ^ 0.44). 

As already mentioned and used for the calibration, the vertex charge Qb can be used 
as fiavour tag for hemispheres where the b- or 6-quark fragmented into a charged b- 
hadron and is thus sensitive to the forward-backward asymmetry A[^. The differential 
asymmetry, computed from the vertex charges in the forward and backward hemispheres, 
versus the direction of the thrust axis is shown in Figure 6, showing good agreement with 
the expectation from the measured value of the pole asymmetry — 0.0982 ± 0.0017 
(from [1]). 

To estimate systematic errors, the following parameters in the simulation and cuts 
have been varied: 

• Lifetimes of 6-hadrons. 

• The oscillation frequency and thus the mixing probability Xd for ^''-mesons. Mixing 
of S°-mesons leaves the contribution from fragmentation particles unchanged but 
reverses the flavour of the weakly-decaying S-meson. The contributions to the vertex 



9 



; 0.14 
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Figure 6: The dependence of the measured differential asymmetry A^-^j on cos 9 thrust for 
the data. A^^jjr has been computed from the vertex charges in the forward and backward 
hemispheres as: A^/^j = (< Q^'^ > - < Qj^^' >)/ < >, where < > is the 
mean of the vertex charge for hemispheres containing a 6-quark. The errors indicated 
are statistical only. The expectation for the pole asymmetry A^^ — 0.0982 (from [1]) is 
superimposed as solid line. 



charge from B - and B -mesons are slightly different. One reason is that i? -mesons 
produce dominantly £)~-mesons whereas B -mesons mainly give L>+, leading to 
different charges at the tertiary charm vertex. 

• The rates of different 6-hadron species, because different 6-hadrons contribute to a 
certain charge (e.g. 5°, B^, A^, S° to Q = 0) and their distributions look slightly 
different. The value of fBu{= /sd) has been set to the measured value in this 
analysis. The values of /b^ and fb-baryon from [1] have been rescaled accordingly to 

ensure fs, + fb-baryon + fsu + /fid = 1- 

• The rates of excited 6-hadrons, namely orbitally excited S**^ ^-mesons and E^*-*- 
baryons. The expected strong decays of these states produce particles (pions or 
kaons) looking partly like fragmentation particles (coming from the primary vertex) 
and partly hke S-decay particles (having large rapidity). Prom the results in [14,15] 
an average Pb-^ = BR{h Bl*^)/BR{h b{u, d)) = 0.24 ± 0.04 is computed. It 

is assumed that Pb** — Pb**^- For E^*^ production DELPHI gives an upper hmit 
in [15]. The following value is taken as a conservative choice^: P^(*) — 0.10 ± 0.05. 

b 

• The rate of so called 'wrong-sign charm' production at the upper W vertex from [16]. 

• The branching ratios of 'right-sign' decays of B^- and i?°-mesons into D - and 



D~-mesons. Because of t{D~) > t{D ), 6-hadron decays with a D~ in the final 
state have charged particles from the tertiary charm decay which are more displaced 
from the primary vertex than in the case of a D , affecting the vertex variables 

'P (,) is defined in the same way as Pb»* . 
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source 


value and variation 


sn% 






(1.542 ± 0.016) ps 


+0.007 


Tb+ 


(1.674 ± 0.018) ps 


-0.129 




(1.461 ± 0.057) ps 


+0.027 


'^b—baryon 


n oris -1- n nc;i t-io 
l^i.zuo ic u.uoij ps 


+0 010 




0.181 ± 0.004 


+0.038 




(8.5 ± 1.3) % 


+0.217 


fb~baryon 


(9.5 ± 2.0) % 


+0.057 




(1.1 ±0.5) % 


-0.187 




0.24 ±0.04 


+0.287 




0.10 ±0.05 


-0.027 


wrong sign charm rate 


(20.0 ± 3.3) % 


-0.001 


BR{B''^+) ^ :D~^°^ + X) 


text 


-0.183 


b fragmentation function 


text 


0.201 


min. # of acc. charged particles 


4,3 


0.217 


Qb cahbration 


text 


0.677 


non-66 background 


±30% 


+0.113 


total 




0.886 



Table 1: Breakdown of systematic errors on f'^. For the 6-hadron fractions, /s^ and 
fb-baryon, their Correlation has been taken into account, for the other sources of sys- 
tematic uncertainties the errors have been added in quadrature. The signs of the er- 
rors given are for an upwards variation of the corresponding physics parameter. For 

SP(S°(+) ^ D~^°^ + X), an upwards variation means increasing BR(B^ — > D + X) 
and adjusting the other, related branching ratios as explained in the text. 

used as input to the neural network. The following branching ratios have been 
used, being consistent with the measured inclusive production rates of -D - and D~- 
mesons in 5-decays: BR{B^ ^ + X) = 15.6%, BR{B^ ^ if + X) = 65.8%, 
BR{B+ ^ D- + X) = 29.3%, BR{B+ ^ D° + X) = 52.1%. These branching 
ratios have been varied by 5% (absolute). The variation has been performed in 
a correlated way to keep the total rate oi D - and £)~-mesons in S-decays and 
BR{B^/B+ ^^ + X) + BR{B^/B+ ^ + X) constant. 

• The 6-quark fragmentation function. The measured function from DELPHI [17] has 
been used. For the systematic error the full difference to the model implemented in 
the simulation (Peterson with = 0.002326) has been taken. This is a conservative 
choice, since the measurement errors are much smaller than the deviations from the 
fragmentation function in the simulation. 

• The cut on the minimum number of accepted charged particles in the hemisphere 
used for the calculation of has been changed from an even to an odd number. 
This could have a systematic effect because charged (neutral) 6-hadrons decay to an 
odd (even) number of charged particles. 

• The calibration of the vertex charge (described in Section 3). The parameters ai 
and 0,2 have been varied independently in a way which results in a displacement of 
the f^'^ versus curve (Figure 4) in the simulation, corresponding to a shift by 
one standard deviation with respect to the data errors. This procedure has been 
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consistently applied also to the 1995 year simulation, even if no correction has been 
applied in this case, because of the good agreement between data and simulation. 
The error from the parameter giving the largest variation in /+ has been chosen. 
• The non-66 background (mainly cc) has been varied by ±30%. 

If not explicitly stated otherwise, the corresponding numbers have been taken from [1]. 
The breakdown of the systematic errors is shown in Table 1. The result, including 
systematic errors, is: 

/+ = (42.09 ± 0.82(stat.) ± 0.89(syst.))%. (4) 



5 Interpretation of the results 

Weakly-decaying neutral 6-hadrons are B^- and S^-mesons, the A^-baryon and the 
strange 6-baryon S^. The rate of charged 6-hadrons is dominated by the S+-meson, 
strange 6-baryons giving only a minor contribution (S^, Jl^). Formally, one gets: 

f f Bd ~^ f Be f b—baryoni 
f — fSu + fb —bar yon' 

(5) 

In [16], the fraction of S^-baryons has been estimated analysing production rates of 
S-r final states: f^- = BR{b ^ S^) = BR{b E') = (1.1 ± 0.5)%. This fraction is 
subtracted from to get the fraction of i?+-mesons in a sample of weakly-decaying b- 
hadrons produced in the fragmentation of 6-quarks (the production fraction of Q^-baryons 
is expected to be negligible). The result is 

/^^ = (40.99 ± 0.82(stat.) ± l.ll(syst.))%, (6) 

where the error from f^- has been added to the systematic error from /"•" taking into 
account the correlation arising from the fact that the uncertainty of is a source of 
uncertainty for /+ (compare Table 1). 



6 Conclusions 

A precise measurement of the production fractions of weakly-decaying charged and 
neutral 6-hadrons has been presented for the first time. The fraction of 6-quarks fragment- 
ing into positively charged weakly-decaying 6-hadrons, and thus the fraction of charged 
6-hadrons in a sample of weakly-decaying 6-hadrons produced in — >■ 66 decays, has been 
measured to be /+ = (42.09 ± 0.82(stat.) ±0.89(syst.))% = (42.09± 1.21)%. Subtracting 
the eI rate (assuming that the rate is neghgible) gives /b„ = (40.99 ± 0.82(stat.) ± 
l.ll(syst.))% = (40.99 ± 1.38)%. This is, so far, the most precise dedicated measurement 
of a production fraction of a specific 6-hadron. The accuracy on fs^ is comparable to the 
accuracy achieved by combining all other information available on 6-hadron production 
fractions. This measurement thus represents significant new information that will impact 
on combinations of measurements to estimate 6-hadron production fractions in jets. 
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